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We give algorithms to factorize large integers in the duality computer. We provide three duality 
algorithms for factorization based on a naive factorization method, the Shor algorithm in quantum 
computing, and the Fermat's method in classical computing. All these algorithms are polynomial 
in the input size. 

PACS numbers: 03.67.Hk,03.65.Ud,03.67.Dd,03.67.-a 



Factorizing large integers has been a great problem in mathematics and has important applications in information 
security. Many efforts have been made to find polynomial algorithms. The most naive method is that we divide 
the input state N by each integer ranging from 2 to \/N, so it needs y/N operations which needs 2^ lo&N steps, an 
exponential function of the input size log N . Until now, there has been no known efficient polynomial algorithm for 
factoring a large integer in classical computing. The most efficient factoring algorithms in classical computing are the 
Fermat's factoring algorithm^. 

With the advent of quantum computer 0, 0, Shor gave a remarkable polynomial-time quantum algorithm for the 
factoring problem, and its computational complexity has been proved to be 0((log 7V)2(loglogiV)(logloglog^V))Q. 
Recently, a new computing model, the duality computer(DC), was proposed based on the general principle of quantum 
interference 5]. The DC exploits the particle wave duality properties of the quantum world, and provides polynomial 
algorithms for NP-complete problems. Finding a marked item from an unsorted database requires only a single step 
in the duality computer, which achieves the holy grail in the unsorted database problem. Recently, the mathematical 
theory of the duality computer has been worked out 0,0. In the DC, there are quantum wave divider and combiner 
operations in addition to the usual unitary quantum gate operations. The divider and combiner operators make 
the connection between duality computing and classical computing very close, and it is possible that many classical 
algorithms can be adopted to the duality computing with exponential speedup. In this paper, we explore this 
connection with the problem of prime factorization as the subject. 

This paper is organized as follows. In Sec^ we briefly review the DC and the implementation of the naive factor- 
ization algorithm. In Sec lIIII we briefly review Shor's algorithm and its implementation in the DC. Compared to the 
quantum computer, the algorithm in the DC acquires some improvement though both algorithms are polynomial ones. 
In Sec lIVI we design two algorithms using DC based on both the Fermat's methods and analyze their complexities. 
The results are discussed and summarized in SecfVl 



The general principle of quantum interference for the quantum system was proposed and based on it, the DC was 
proposed The DC makes use of not only the particle nature but also the wave nature of microscopic objects. The 
quantum wave of the DC are split into several sub- waves by the quantum wave divider(QWD) and they pass through 
different routes where different computing gate operations are performed. These sub-waves are then re-combined to 
interfere to give the computational results at the quantum wave combiner (QWC). 

The naive factorization algorithm can be easily implemented using the Long algorithm for unsorted database search 
(BJ- We need a two register DC. The DC algorithm is given as follows 

Step 0: foundlist=[null] 

Step 1: Prepare the the state of the first DC register in the equally distributed state and the second register in the 
state, 



I. INTRODUCTION 



II. THE NAIVE FACTORIZATION DC ALGORITHM 




(1) 



where [k] means the integer nearest to k. 

Step 2: Let the DC go through QWD, so that it divides the wave into two sub-waves, 



Wu) = 2 y ~gw\ (|2> + 13> + ' ' ' + lPi) + " ' + 1[ 1 + 1>)|0> ' (2) 



\<p d ) = ^^(12) + |3> + ■ ■ ■ + \Pi) +-- + \[VN\ + i»io>. 
Step 3: Perform the following function to the lower-path sub-wave, 



, 1 if N mod i = , , 

W*) = 1 otherwise (3) 



Then the state of the lower-path becomes 



y d ) = I^_L(|2>|0) + | 3 )|0) + ■ ■ ■ + |Pi>|l) + • • • + \Wn\ + i»|o». 

Step 4: Apply a query to the lower-path wave: if the second register of the item is 1 retain the sign of the item, 
reverse the sign of other items and the items in the foundlist. 

Step 5: Recombine the two sub-waves \ip u ) and \ip' d ) at the QWC. Then \ipf) = \pi) + \p 2 ) + • • • + \p m ) apart from 
a normalization constant, namely \ipf) contains all the factors except those in the foundlist. 

Step 6: Make a read-out measurement, and one of the unfound prime factors p» is found out with probability 1. If 
nothing is found, then go to step 8. 

Step 7: Add pi to foundlist. 

Step 8: Stop. 

If N = pq, then the above algorithm runs from step to step 7 to give the smaller factor, and then it goes to step 
and flows to step 6, and then goes to step 8 and stops. 



III. THE DC FACTORIZATION ALGORITHM BASED ON THE SHOR ALGORITHM 

The Shor algorithm finds the period of the function f a _N = a x mod N. In the quantum computer, the registers of 
the quantum computer is divided into two, and the Shor algorithm contains the following steps: 
Step 1: Prepare the first register in the evenly distributed state and the second register in state 



IV'i) = ^I>>|0>, ( 4 ) 



o 



where N 2 < q < 2N 2 . 

Step 2: Apply the function f a jf to the state and store the results in the second register, 



IV2) = \ ~J2\i)\a x mod N). (5) 
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Step 3: Apply a Fourier transform to the first register. 

Step 4: Makes a measurement on the first register and an integer will be found. 

Repeat steps 1 to 4 a sufficient number of times and from the measured results a period is inferred, from which the 
period of the function f a ,N(%) is determined. 

In the DC there is no need to apply the Fourier transform on the first register. Here we present a modified Shor 
algorithm in the DC with N = 21 and a = 2 as an example. Taking q = 512, after applying the function / 0) jv on the 
DC, we let wave goes through the QWD, so that its wave function is split into two parts, and in the lower path we 



reverse the sign of all items except those whose second register number is 1 and the first register is not zero, namely 
the upper and lower path wave functions are now 

\tp u ) = 5^=(|0>|1) + |1)|2) + |2)|4> + |3)|8) + |4)|16) + |5)|11> + 

|6)|1) + |7)|2) + |8)|4) + |9)|8) + |10)|16) + |11)|11)+ (6) 
|12)|1) + ---), 

\ Vd ) = \-±={ - |0)|1) - |1)|2) - |2)|4) - |3>|8> - |4)|16) - |5>|11> + 
* V 512 

|6)|1) - |7)|2) - |8>|4> - |9)|8) - |10>|16> - (?) 
|12>|1> ), 

respectively. 

Next step, we combine the sub-waves at the QWC, and apart from an normalization constant the wave becomes 

\<p f ) = (|6)|l) + |12)|l) + ... + |504)|l) + |510)|l)) 

= (|6) + |12) + ... + |504) + |510))|l). (8) 

Make a read-out measurement of the first register, one will obtain one of the following numbers: 6, 12, 18, 24 ... 
504, 510 with equal probability. Repeat the procedure a number of times, the smallest increment will be found and 
hence the required period. 

IV. DC FACTORIZATION ALGORITHM BASED ON FERMAT'S METHOD 

On classical computers, Fermat's methodQ] for factoring an odd integer N — pq where q, and p are prime integers 
consists of finding positive integers X and Y that satisfies N = x 2 — y 2 . The key points lie in finding those X and 
Y satisfying N = X 2 — Y 2 so that we can get the result p = (X + Y)/2,q = (X — Y)/2. Fermat's method is quite 
efficient if p/q is near 1, but less efficient if p/q is not near 1. The computational complexity for Fermat's method is 
0(q - \/N), and in worst conditions is 0(N/2 - y/W) 

We design an algorithm under worst conditions for a DC based on Fermat's method and analyze its complexity. 
We denote the integer bigger than or equal to (N/2 — \/N) as M 

1. Make the initial state 

1 M 1 1 

l*°> = 77^ E l*> = rm^ + 1* 1 * + • • • + 1*")) = + ]VN+ 1} + ' ' ' + lN/2)) (9) 

i=0 



2. Let \tpo) go through the QWD and it is divided into: 

N/2 



lv-> = irm E \**> = + ^ + ' ' ' + 1^)) = ^ + 1^ + !> + • • • + \n/2)) (io) 

»=o 



M 



\<P*) = 77T7 E 1^) = irm^) + I* 1 ) + ' ■ ■ + 1^)) = TTfjd^O + 1^ + 1) + ■ ■ • + W/2)) (11) 
ly M ~^ 2v M zv M 

3. Apply on the upper-sub wave the following query function f Xt N 



f _ J +1 if \/ x 2 — N is an integer , 10 >, 
ix ' N ~ \ -1 otherwise (12j 

4. Let the two sub-waves go through the QWC and it is combined into \X), then we can make a measurement to 
read out this X, and the Y is then y/X 2 — N. Then p and q are obtained. 

The complexity for our algorithm is C(log (N/2)) 



V. SUMMARY 



In this paper, we propose three duality algorithms to factorize large integer number. It has been shown that Shor's 
algorithm in a duality computer can be simplified a little. The naive factorization algorithm can be implemented in 
a duality computer polynomially. We also show that the Fermat's method in classical computing can be adopted in 
duality computing with a computational complexity of O(logiV). 

Compared to quantum computing, the computational complexity of prime-factorization algorithm is the same which 
is not remarkable. The important point shown in this work is that exponential algorithms in classical computing can 
be implemented in duality computer with a polynomial computational complexity. This provides a convenient avenue 
in constructing duality computing algorithms. 
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